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Abstract 

In this paper we study a queue with Levy input, without imposing any a priori assumption on 
the jumps being one-sided. The focus is on computing the transforms of all sorts of quantities 
related to the transient workload, assuming the workload is in stationarity at time 0. The results 
are simple expressions that are in terms of the bivariate Laplace exponents of ladder processes. 
In particular, we derive the transform of the minimum workload attained over an exponentially 
distributed interval. 
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1 Introduction 

In this short communication we focus on the analysis of various random quantities related to the 
transient of a Levy-driven queue. Throughout, we denote the driving Levy process by Xt, t > 0, 
uniquely characterized through its Levy exponent 

= LogEe'^^i, 6 »gM. 

As usual, the corresponding queueing proeess (or: workload process) Qt, t > 0 is defined as the solution 
to the Skorokhod problem [5, Ch. II], i.e., 

Qt = Qo + Xt — (A^ -|- Qo) A 0 = Qo V (—A^) + Xt, t > 0, (I) 

where A^ := mf^gjo^tj A^ is the running minimum and Qo is assumed to be independent of the evolution 
of the driving process Xt for t > 0; similarly we let At := sup^g^ tj A^. 

It is well known that Qt converges to a stationary distribution as t —)■ oo if and only if At has a negative 
drift, that is, limt-too Jft = —co a.s.; we therefore assume that this condition is in place throughout 
this work. If EAi is well-defined, then this condition is equivalent to requiring that EAi G [—oo,0), 
see [6, Thm. 7.2]. It is also a standard result that the stationary distribution of Qt coincides with the 


I 



distribution of X := Xoo, the overall supremum; this is a property that is often attributed to Reich 
[5, Eqn. (2.5)]. 

Importantly, in the sequel we systematically assume that Qq has the stationary distribution] as will 
become clear, this assumption plays a crucial role in the analysis. As mentioned above, the primary 
objective of this note is to study various transient metrics related to the process Qt- All our results are 
in terms of transforms of the quantities of interest, in addition transformed with respect to time. In this 
respect, it is recalled that taking transforms with respect to time essentially amounts to considering 
t = eg, where eg is an exponential variable with rate q (i.e., with mean q~^), sampled independently 
of everything else. 

This note is organized as follows. In Section 2 we sketch preliminaries, related to splitting at extrema 
and the Wiener-Hopf factorization. In Section 3 we focus on the distribution of the minimum workload 
attained over an exponentially distributed amount of time. Section 4 provides transforms of some 
other quantities of interest, distinguishing between two cases: (i) the ongoing busy period and (ii) the 
finished initial busy period, in which case the focus is on the so-called unused capacity. Finally, 
Section 5 focuses on the minimal workload in a queue conditioned to be positive. We conclude the 
paper by mentioning some challenging open problems. 


2 Preliminaries 


Let us state some well-known facts about splitting and the Wiener-Hopf factorization, all of which 
can be found in e.g. [1, Ch. VI]. Consider the process X on the interval [0, eg] and let 

Ge^ := sup{t G [0, eg] : Xt A Xt- = 

be the (last) time of the minimum. It is well known that the process X splits at its minimum, see [I, 
Lem. VL6]. This result can be stated in a simple form when considering the process X[ := Xt for all 
t / and X'q ;= Vg_^. Namely, the processes 

{X[, t G [0, GgJ} and - X^^^ , t G [0, e, - GgJ} 

are independent (note that X' and X differ only if X jumps up at Gg — this jump should be included 
in the right process, but not in the left). In particular, the above implies that X^^ and Ag^ — Ag^ are 
independent. Splitting in the queueing context used earlier; see e.g. [3]. 

Let us define the following two functions for a,/I G C with il?(a), il?(/3) > 0; 

K{a,(3) = exp if f (e“* — e“"*“^*)-P(At G dx)dt ) , 

\Jo 2(0,oo) i ) 

n{a,j3) = exp if f (e“* — e“"*“^^)-P(At G dx)dt ) , (2) 

\Jo J{-oo,0] t J 

which are the Laplace exponents of the so-called (strictly) ascending and (weakly) descending bivariate 
ladder processes, respectively. In general these functions are defined up to a multiplicative constant; 
the trivial scaling chosen in (2) implies that 


K{a,0)K{a,0) = exp 



= a, 


a > 0, 


(3) 
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by virtue of the Frullani integral identity. This particular scaling simplifies the calculations in the 
proofs, but it does not affect the results of the paper. The Laplace exponents can then be used to 
express the so-called Wiener-Hopf factors: 


K{q + «,/?)’ 

K{q + a, /?) ’ 


where := inf{t G [0,6^] : Xt V Xt- = Xg^} is the (first) time of the maximum. 
Finally, for /3 G M we have 


q 

q-fjiP) 




K{q,0) K{q,0) 

K{q,il3) K{q, -i/3) ’ 


which according to (3) yields 


q - = K{q^ i/3)K{q, -i/3). 


(4) 

(5) 


( 6 ) 


Let us finally remark that one has to distinguish between first and last extrema only in the case when 
X is a compound Poisson process; the same is true also for open and closed integration sets in (2). 


3 Minimum workload 

In [4] the distribution of the minimum workload during [0, e^) was characterized for the case the 
driving Levy process is spectrally one-sided; it was assumed that the workload is in stationarity at 
time 0. In this section, this result is extended to the spectrally two-sided case. Several ramifications 
of this result are presented as well. 

Throughout we let be the running minimum of the workload, i.e., is the minimum of Qs over 
s G [0,t], and we assume that Qo has the stationary workload distribution. 


3.1 Transform of minimum workload 

This subsection shows how to express the transform of in terms of the Wiener-Hopf factors. Observe 
that 

Q, = (Qo + M) V 0 = X, + Qo V i-x,) (7) 

which follows by considering the event Qo > —2Lt its complement separately. To this end, we first 
rewrite the definition of Q*, as given in (1), in terms of Xt — X_t and Q^: 

Qt = {Xt - Xt) + X, + Qo V i-Xt) = {Xt - Xt) + Q^, 

which is also easily seen from Figure 1. This means that we in particular have 


Qe, = (Xe,-X,J + Q^^, 

where the two terms on the right are independent, because of (7) and splitting at the infimum, i.e. 
Xg^ — and Xg_^ are independent, see Section 2. Using (5) we thus obtain the identity 




K{q,e) 


( 8 ) 


for 0 > 0. Since Qo has the stationary workload distribution, so does Qg^. As a consequence, Qe, has 
the distribution of X, which yields the following result, see also (5). 
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Proposition 1. For any 6,q>0 


7c(0, 9) K{q, 0) 


(9) 


It is noted that Equation (9) generalizes the Endings of [4] which just cover the spectrally one-sided 
cases; it is readily verified that plugging in expressions for F for the one-sided cases the formulas 
obtained are in accordance with those derived in [4], 


Remark 1. The transform of Q has a simple explicit form also in the case when Qo is exponential 

-EfJ 

(with mean X~^): 

E e~^^o = P < eg) = P (Qq + = P(Qo < eg - = 1 - E 

^ ^ ^ g(g,0) 

X + 6 K{q, X) 


and then (8) yields Ee . Throughout this work we assume that Qo has the stationary distribution, 
but virtually all results carry over to the case of the exponential initial distribution. 


We conclude this section by deriving a related, elegant identity. To this end, we dehne by 


r := inf{t > 0 : Qt A Qt- = 0} = inf{t > 0 : Qt = 0} a.s.. 


the time when the system becomes empty for the first time; r is also referred to as the residual busy 
period as seen from time 0, and the equality follows from [1, Prop. VI.4]. Clearly, P(r = e^) = 0 and 
as a consequence 

P = o) = P(r < eQ. 

Taking 0 —)■ oo in (9), and noting that K{q,9)/K{0,6) —)■ I which follows from (2), we obtain 

Ee-^^ = P(r < eg) = P = Oj = = Ee"''^ (10) 

showing that r and G := Goo, the (first) time of the overall maxima, have the same distribution (which 
was concluded, using another line of argumentation, in [3] as well). 


4 More refined quantities 

In this section we analyze a few more refined quantities that are related to the minimum workload. 
In this respect observe that 

Ge, = sup G [0, eg] : Qt A Qt- = 

is also the (last) time of the minimal workload, see Figure I. In the sequel we find it convenient to 
distinguish between the following two cases: r > Cg and r < Cg, i.e., if the residual busy period is 
still going on at time Cg, or not. In the second case we look at the first and the last times when the 
workload attains the value 0, i.e. the time r when the hrst busy period finishes and the time G^,^ when 
the last busy period starts. 
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Figure 1; Schematic queueing process 


4.1 Ongoing busy period: t > Cq 

Note that (9) combined with (10) provides us with the equality 

E ^—‘^ 1 ; T > Cq'j = Ee — P(r < Cq) 

^ 7t(0, 0) K{q, 9) k{0, 0) _ 7t(0, 0)(K(g, 6) - 7t(0, 9)) 

7t(0,0) 0) li{q,0) K{0,9)K{q,0) ’ 

describing the distribution of Q on the event t > Cq. The following result gives a simple expression 

^<2 

of the joint transform of various quantities related to the minimal workload in an ongoing busy period. 


Proposition 2. For any 6,a,j3,'y,q > 0, 


E e“ -'=9^ « 


;T> Cq] = 


q K{0,0){K{q + f3,6 + a) — k{ 0,9 + a)) 


P + q k{Q, 9 + a)K{q + ^,a) 

Proof. By appealing to the usual splitting argument, we can write 


E e -"9 


T > e„ ) = E ( e ; T > e„ ) Ee 


-Ke,-n) 


^ / -SQ 
= E e —'=9; T > e. 


K{q, 0) k( 0, 0)(K(q, 9) - k(0, 9)) 


</ T. 


see (5) and (11). Now realize that 


E f(ea)e-P^‘^ = 


-9eA _ _J_ 


P + q 


K{q + /3,0) K{0,9)K{q + /3,0) 


E/K+/3) (12) 


for any Borel function /. It therefore immediately follows that 




E (^e ; r > Cg ) = 

Replacing consistently q + P by q, and j3 by —/3, we thus obtain 


\ Q ir-i f ~^Q J^g \ 


E 



T > Cc 


q K{0,0){K{q +13,9) - k{ 0,9)) 
I3 + q K{0,9)'K{q,0) 
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Finally, again relying on the splitting argument, we arrive at 




= E e 


-{a+e)Q -PG^ 




q k(0 , 0)(K(g +/?, 0 + a) — k(0 , 0 + a)) k((?, 0 ) 


P + q 

see (5), completing the proof. 


«;(0, 6 + a)K{q, 0) 


K(g + 7,a)’ 


□ 


4.2 Finished busy period: t < Cq 

In the case t < eg we define two random quantities: 


Ue,:=-{Qo+Ke,), D:=-{Qo + Xr), 

see Figure 1(b). The first can be interpreted as the unused capacity, and the second as the sudden 
unused capacity at the end of the initial busy period. The latter quantity may be of limited interest, 
but it can be included in the joint transform without any additional work. First we start with a basic 
result complementing (11). 


Proposition 3. For any 6,q>t), 

E ; T < e,) = na.O)^) - 

V ) K(g,0) n{q,0) 

Proof. First we compute for 9 = iR 

E ; r < eg) = - E ; r > e^) 

k( 0,0) /i((7,0) K{0,0){K{q,—6) — k{0, —9)) 

K{0,—9)K{q,9) k{ 0,—9)F{q,0) 

see (5), (4) and (11). Using (3) this can be rewritten as 


«^( 0 , 0 ) f q 

K{q,0)K{0,-9) \Kiq,9) 


K-{q, -0) + «^( 0 , -9) 


Finally, using (6) we can express K{q,—9) = {q — 'ip{—i9))/K{0,9). Plugging this in yields that the 
expression in the previous display equals 

75(0, 0)k(0, 0) f Tp{—i9) 'ip{—i9)\ 

-F{q, 0)V’(-i6l) V K{q, 0) k(0, 9) ) ’ 

which easily reduces to the expression in the statement. Finally, analytic continuation shows that it 
is true for all 6* G C with iR.{9) >0. □ 


Proposition 4. For any 9, a, j3, 7 , u,v,q > 0, 

q k( 0 , 0 )(K(g + tt, a +/3) — 7$(0, a +/?)) 
q + u K{q + w,j)K{q + v,(3) 
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Proof. By splitting we have 


E < e,) = EV"<'--a,) = «( 0 . 0 )fefa.«)-il(°.«)) 

V Hj \ H/ k(( 7 , q;)k ((7 + n, 0) 

see also (5) and Proposition 3. Moreover, using (12) we write 




As before, changing parameters leads to 


E = —-—E "^'=9+i*+“(®9+“ r < Cg+^t^ 


V ' V g + tt V 

q 75(0,0)(K(g + u, a) — k( 0, a)) 

g + tt K(g + tt, a)75((7, 0) 

Next, we apply the strong Markov property at the time r to obtain the identity 


^ / aUe, ; r < eg ) = E U-^D-ur 


^ = E (e-“^-“^; r < eg) E , 


which leads to 


—OiD—ur 


; 'T < eg) = 


/t(0, 0 )(k ((7 + u,a) — k{0, a)) 
q + u 


using (4) and (2). Finally, we write using the strong Markov property at r and splitting at Gg 


-aD-l3Ue„-'yQe„-uT-v{G -T)-w{eq-G^ ) 


■JU T < e. 


= E < eg^ Ee^-"9 ’'-'=^Ee 

F(0, 0)(K(g + u, a +/?) — k(0 , a +/?)) K{q,0) K{q,0) 

q + u K{q + v,(3)K{q + w,'y)' 

which in view of (3) completes the proof. 


5 Minimal workload in a queue conditioned to stay positive 

In this section we focus on the law of Q conditional on the queue not having idled between 0 and 
eg, i.e., r > eg, and provide an alternative representation of this law in the limiting case when ( 7 1 0- 
We also comment on the relation of this limit law to Levy processes conditioned to stay positive. 

It follows directly from (10) and (11) that 

E 1. > . E . > ,)/p,. > (n, 

Note, however, that the corresponding conditional law does not have a direct link (via Laplace trans¬ 
form) to its transient counterpart, i.e. when Cg is replaced by t. 

In the following we assume that 75'(0, 0) < 00 (and so also 75'(0, 6) < 00 for 6* > 0), where the derivative 
of 75(g, 6) is taken with respect to q. According to (10) this requirement is equivalent to Er = EG < 00 ; 
see Proposition 6 for an example when this assumption does not hold. Now it follows from (13) that 

limE (^“3-. It > e,) ^ ^ !°"!!!°'"!!) (14) 

(?to V k(0, 0) k'( 0, 0) log(«;(0,0))' 

Along the same lines, one can establish the generalization 

/ I ^ ^ = 75(0,0) log(75(0,6* + Q;))' 

V 75(0,0;) log(75(0,0))' 


limE 

qiO 
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Remark 2. It must be noted that the limit law of Q \t > Cq as q f 0 does not eoincide with the 

— 

so-ealled quasi-stationary distribution of the minimal workload, i.e. with that of Q^\t > t as t ^ oo, 
even though —)■ oo a.s. as g | 0. Roughly speaking, conditioning on {r > eg} does not only force r 
to be large, but also makes Cg to appear smaller. In this respeet it may be helpful to mention that 


lim E ( e I r > e, 
qJ-O V 


1 -Ee-®^ 


i.e. the limit law of eq\T > eg is a proper distribution, the residual life distribution assoeiated to t. For 
related results on the quasi-stationary behaviour of refleeted one-sided Levy processes we refer to e.g. 
f7j. 


Remark 3. Importantly, one way to construct a Levy process started in x and eonditioned to stay 

positive (in the usual sense) is to eondition on {r > Cg} and then let (7 0, see [2, Prop. 1]. The 

distribution of the infimum of this process is charaeterized in [2, Thm. 1]. Note, however, that the 

limit distribution of Q |r > e„ can not be obtained from this result by integrating with respect to 
- 

E((5o £ dx). The main reason is that 


, -aQ , , °'—‘^i:T>eq) , EEo„(e ®9 ;r>e„) -aQ , , 

EEQ„(e ^^flT>eg)=E ^ = E(e > Cg), 


^Qoi'^ > eg) 


EPqo(^ > eg) 


beeause the event upon whieh we condition depends on Qq. 


Proposition 5. Assume that k'(0,0) < 00 . Then the limit laws of 


Q \t > Cg, and Xg |Xe„ > 0 

eoincide as g | 0 . 

Proof. Using (2) observe that 

q log{K{q, ~ J j qe~‘^^~^^E {Xt G dx) dt = E ; Xg^ > 0^ 


Hence 


= lim = lim ^ E (e-*-''-. | JC, > o) , 

log(/t(0,0))' gto log(K(g, 0))' gto P(Xg > 0 ) gto V / 


log(/t(0,0))' gto log(K(g, 0))' gto P(Xg^>0) 
The proof is complete in view of (14). 


□ 


This result has a simple intuitive explanation. Firstly, on the left hand side we have the limit of 
Qo + XgJQo + Xg^ > 0. Secondly, it holds that Xg^ = (Xg^ — Xg^) + Xg^, where the two terms on 
the right are independent and the distribution of the first converges to that of Qo as g | 0 . 

Let us conclude by giving yet another representation of the limit law of Q |r > e„ in the spectrally 

—e? 

one-sided cases. 


Proposition 6. Assume that X is either spectrally positive or speetrally negative. Then the limiting 
law of Q \t > Cg as q f H is the residual life distribution associated to Qq, i.e. 

-Eg 


/ _ 

lim E I e I r > Cg 
gto V ^ 


1 — E e 

0 EQo 


(15) 


Moreover, k'(0,0) = 00 if and only if X is a speetrally positive process with var(Xi) = 00 , in which 
case EQo = co and Q \t > Cg converges to 00 as g | 0. 

- Gq ^ 



Proof. The identity (15) is easily verified using (14) and the explicit expressions for k(-, •) in both 
cases, see, e.g., [6, Sec. 6.5.2], These explicit expressions also show that in the spectrally negative case 
we must have 7«'(0, 0) < oo, whereas in the spectrally positive case 7«'(0, 0) < oo iff 4>”{0) = oo, where 
4>{9) = \og'Ee~^^^. The latter is equivalent to var(Xi) = oo and implies that (15) results in 0. □ 

It seems unlikely that (15) holds in general. It would be interesting to characterize all the Levy 
processes drifting to —oo for which (15) is true. Another challenging problem is to express the joint 
transform (or, closely related, the joint transform Ee “'*«) through the functions 

k{-, •) and k{-, •). 
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